SUBROUTINE FFT2D(A,N,M) 
    DIMENSION A(N,M)
    COMMON/WSPACE/WORK(2048)
    CALL FFTINIT(N,WORK)
    DO 100 J=1,N
100 CALL FFTCOLL(A(I,J),N,1,WORK)
    CALL FFTINIT(M,WORK)
    DO 200 I=1,N
200 CALL FFTROW(A(I,1),M,N,WORK)
    RETURN
    END